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Context. The B/C ratio is used in cosmic- ray physics to constrain the transport parameters. However, from the same 
PjQ . set of data, the various published values show a puzzling large scatter of these parameters. 

Aims. We investigate the impact of using different inputs (gas density and hydrogen fraction in the Galactic disc, 
source spectral shape, low-energy dependence of the diffusion coefficient, and nuclear fragmentation cross-sections) on 
the best-fit values of the transport parameters. We quantify the systematics produced when varying these inputs, and 
O |i compare them to statistical uncertainties. We discuss the consequences for the slope of the diffusion coefficient 5. 

I ■ Methods. The analysis relies on the propagation code USINE interfaced with the Minuit minimisation routines. 

Results. We find the typical systematic uncertainties to be greater than the statistical ones. The several published 
values of 5 (from 0.3 to 0.8) can be recovered when varying the low-energy shape of the diffusion coefficient and the 
convective wind strength. Models including a convective wind are characterised by S > 0.6, which cannot be reconciled 
with the expected theoretical values (1/3 and 1/2). However, from a statistical point of view (x 2 analysis), models 
with both reacceleration and convection — hence large S — are favoured. The next favoured models in line yield S, which 
can be accommodated with 1/3 and 1/2, but require a strong upturn of the diffusion coefficient at low energy (and no 
convection) . 

Conclusions. To date, using the best statistical tools, the transport parameter determination is still plagued by many 
unknowns at low energy (~ GeV/n). To disentangle all these configurations, measurements of the B/C ratio at TeV/n 
energies and/or combination with other secondary-to- primary ratios is necessary. 
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. 1. Introduction ~ 10 — 20%. Although such precision is expected given the 

,— ' ' . _ . , nrlT ,, . accuracy of the data, one may wonder whether these re- 

• ■ A central question in Galactic cosmic ray (GCR) physics ^ are nQt dominated b syste matic uncertainties from 

£ . J*> m tr f determination of the transport parameters. ^ j in dients and assumptions made to do the cal- 

. I he standard procedure consists m fatting a secondary-to- cu l a tion 
;_i ■ primary ratio, e.g. B/C. In general, the parameters derived 

C3 in different studies provide inconsistent values, especially This is an important issue since in our pre- 
fer the slope 8 of the diffusion coefficient. It is important vious studies, the best-fit slope observed in a 
to understand the origin of these differences. diffusion/ constant-convection/reacceleration propagation 
The first attempts to get the best-fit value of the trans- model (|Maurin et al.||200lj |200J)— confirmed by our recent 
port parameters as well as th e ir sta tistica l uncertainties MCMC analysis (Paper II)— points to 8 w 0.75 - 0.85. 
were carried in iMaurin et al.l (12001 and iLionetto et all This is at variance with the result of the GALPROP 
((20051) . Recently, we implemented a Markov chain Monte code ([Strong fc Moskalenkoj [1998|), where the best-fit 
Carlo (MCMC) to estimate the probability-density function models correspond to smaller 8 « 0.3 - 0.4 (|Lionetto et al.| 
(PDF) of th e transport and source p arameters for Galactic l2QQ5|) in either a diffusion/lincar-convection or a diffu- 
cosmic rays (|Putze et all20Ml20ia hereafter Papers I and sion/reacceleration model. Moreover, none of these results 
II) . This sound statistical technique applied to current B /C a g ree with the best-fit slopes from leaky-box models (with 
data allowed us to estimate the statistical uncertainties rigidity cut-off), where 8 ra 0.5 - 0.6 (e.g., |Webber et al.| 
on various parameters. We found typical uncertainties of 12Q0J— no reacceleration—, or Paper I— with reaccelera- 

tion). Such a scatter in 8 was also observed in IJones et al.l 

Send offprint requests to: D. Maurin, dmaurin@lpnhe.in2p3.fr (|200ll ). 
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The propagation code USINE (Maurin et al., in prepa- 
ration) provides the GCR fluxes in both the framework 
of the leaky-box model (LBM) and the diffusion model 
(DM). Associated with an efficient minimisation tool for 
finding the best- fit parameters of a model (w/wo convec- 
tion, w/wo diffusive reacceleration) , it allows these differ- 
ences to be addressed thoroughly. We also investigate how 
sensitive these best-fit parameters are to various input in- 
gredients/parameters. 

The paper is organised as follows. In Sect. [2J the DM 
used is briefly recalled. In Sect. [3J the free parameters of 
the study, our methodology, and the input ingredients — the 
systematic effects of which are studied in this paper — are 
described. In Sect.|H we discuss the consequences of varying 
the gas characteristics on the transport parameter determi- 
nation. In Sect. [5J we focus on the effect of the source spec- 
trum on the best-fit values for S in models with or without 
convection/reacceleration. In Sect. [SJ a similar analysis is 
carried out for the low-energy shape of the diffusion coeffi- 
cient. We repeat again the analysis in Sect. using various 
production cross-section sets. In a final step, we explore in 
Sect.[8] the effect of biasing B/C HEAO-3 data at high en- 
ergy. We summarise, discuss our results, and conclude in 
Sect. ED 

2. Description of the diffusion model 

The models and the equations are described in Paper II, to 
which we refer the reader for a complete discussion. 

2.1. Diffusion equation 

The differential density of the nucleus j is a function 
of the total energy E and the position r in the Galaxy. 
Assuming steady-state, the transport equation can be writ- 
ten in a compact form as 

d ( ■ ■ dN j \ 
£^ + -(^-c>— J=^. (1) 

The operator C (we omit the superscript j) describes the 
diffusion K (r, E) and convection V(r) in the Galaxy, the 
decay rate T ia d(E) = 1/(7x0) for radioactive species, and 
the destruction rate T ine i(r,E) = J2ism n iSM(r)vo- mc i(E) 
on the interstellar matter (ISM): 

C(r, E) = -V- (KV) + V • v + r rad + r incl . (2) 

The coefficients b and c are respectively first and second 
order gains/losses in energy, with 

, / ™ / dE 1 V.F /2m + E k \ 

b ^ = <lt>^-— E *(lZ+l£) (3) 

+ {1 + P2) xK 

E p ' 

c(r,E) = /3 2 xK pp . (4) 

The coefficient K pp is the diffusion coefficient in momentum 
space, and it can take several forms (see later). 

2.2. Geometry of the Galaxy 

The Galaxy is modelled to be a thin disc of half-thickness 
h, which contains the gas and the sources of CRs. This disc 



is embedded in a cylindrical diffusive halo of half-thickness 
L (where the gas density is assumed to be equal to 0). A 
constant wind V(r) = sign(z)-V c xe z , perpendicular to the 
Galactic plane, is assumed. In this framework, CRs diffuse 
in the disc and in the halo independently of their position. 
These assumptions allow for semi-analytical solutions of the 
transport equation, as the interactions (destruction, spalla- 
tions, energy gain and losses) are restricted to the thin disc. 
Such semi-analytical models reprod uce all salient features 
of fu ll numerical approaches (e.g., IStrong fc Moskalenkol 
fl998h . 

In this study, the disc half- height is set to h = 100 pc. It 
is not a physical parameter per se in the thin-disc approx- 
imation, but the phenomena occurring in the thin disc are 
related to it . The physical parameter is the surface density 
S: should a different h value be used, a rescaling always 
allows obtaining the same E. 

Considering the radial extension R of the Galaxy to be 
infinite leads to the ID ve rsion of the DM. This geometry, 
used in I Jones et al.l (|2001[) and in our Paper II, is also used 
in this analysis. The corresponding sets of equations and 
their solutions are presented in the Appendix of Paper II. 
They are not repeated here. 

3. Methodology 

As in Papers I and II, three different classes of diffusion 
models are considered. In addition, for completeness, we 
also treat the pure diffusion case: 

— Model = {K , 5}, i.e. pure diffusion (V a = V c = 0); 

— Model I = {Kq, 8 V c }, i.e. no reacceleration (V a = 0); 

— Model II = {Kq, 5, V a }, i.e. no convection (V c = 0): 

— Model III = {K , 6, V c , V a ). 

The parameters Kq and S come from the standard form 
assumed for the diffusion coefficient, namely, 

K(E)=pK 1l 5 . (5) 
3.1. Constrained parameters 

At most, for Model III, four transport parameters need to 
be determined from CR data: 

— A'o, the normalisation of the diffusion coefficient (in unit 
of kpc 2 Myr- 1 ); 

— S, the slope of the diffusion coefficient; 

— V c , the constant convective wind perpendicular to the 
disc (km s _1 ); 

— V a , the Alfvenic speed (km s _1 ) regulating the reaccel- 
eration strength [see Eq. ((5])] . 

As in other studies, we use the B/C ratio to constrain these 
parameters. In DMs, the halo size of the Galaxy L is an ex- 
tra free parameter. It cannot be determined solely from the 
B/C ratio because of the well-known degeneracy between 
K and L. In this study, we choose to fix L. The value of 
the transport parameters for any other values of L can be 
obtained from simple scaling laws, as presented in Fig. 5 of 
Paper II (note that 5 does not depend on L). To keep the 
discussion as simple as possible, this paper is based on B/C 
data alone, with four free parameters {Kq, 6, V c , V a } and 
the halo size set to L = 4 kpc. 
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Table 1. Reference values for the inputs. 



Input name 



Default value/dependence/set 



Gas 


E = 6.17 ■ 10 20 cm 2 , / H = 


90% 


Source spectrum 


??s = -l, a + 5 = 2.65 




K(E) and K PP {E) 


Slab Alfven (SA): Eqs. © 




Cross-sections 


W03 (Webber et al. 20031 




Data 


B/C, dataset F f 





t 31 data points from IMP7-8, Voyager 1&2, ACE, HEAO-3, 
Spacelab, and CREAM04 

Note 1. In subsequent tables of the paper, all results obtained with 
the default inputs are in italics. 



3.2. Inputs and default configuration 

The inputs that we examine and vary are the following 
(details and references are given in Paper II, Sect. 2.4): 

— gas surface density Sism = 2/misM and hydrogen frac- 
tion /h (in number), 

— source spectrum Q^(E) = qj/3 vs R~ a (qj is scaled to 
match the measured elemental fluxes at ~ 10 GeV/n), 

— spatial/momentum diffusion coefficie nts using different 
turbulence models (|Schlickc iscr 20Q2), 

— production cross-sections. 

The reference (default) values for these inputs are gathered 
in Table [1] The default diffusion coefficient K pp (in mo- 
mentum space) is taken from the mode l of minimal reaccel- 
eration by the interstellar turbulence (jOsborne fc Ptuskinl 
[l98llSeo fc Ptuskinl[l99l : 



K rm x K 



<5 (4 - 5 2 ) (4 — (5) ' 



(6) 



where V a is the Alfvcnic speed in the medium. 

Concerning the experimental data used to fit the B/C 
ratio, as in Paper II, we use the following dataset (de- 
noted F in Paper II), which consists o f 31 B/C data points 
(see Fig. [5] of this pa per): IMP7-8 (iGarcia-Munoz et al.l 
Il987ft. Vovaeer 1&2 (iLukasiak et al.l 11999ft. ACE-CRIS 



( de Nolfo et al.l 12006ft. HEAO-3 (lE ngelmann et~aTl I1990D 



Spac elab (ISwordv et all 11990ft . and CREAM (|Ahn et all 

joog) . 



3.3. Minimisation routine 

In Papers I and II, we adapted the MCMC technique for a 
propagation code. This allows the PDF of the parameters 
to be obtained along with their statistical uncertainties. 
However, this technique relies on thousands of calculations 
for a single setting of the inputs, which is not optimal for 
speed. Here, we are only interested in the best-fit values, not 
in the PDF of the parameters. Therefore, we use the Minuit 
library (a CERN library), which provides minimisation rou- 
tines. Instead of a few hours of distributed calculations for 
the MCMC technique, a few minutes on a workstation are 
enough to obtain the best-fit parameters. 

A few configurations have been cross-checked with the 
MCMC technique, to ensure that the typical widths of the 
PDFs remain the same, whatever the input ingredients used 
in the calculation. This allows us, for a given propagation 



model calculated from different input ingredients, to com- 
pare the resulting scatter in the transport parameter val- 
ues — systematic uncertainties (hereafter SystUnc) — to the 
typical statistical uncertainties (hereafter StatUnc) found 
with the MCMC technique. 



4. Influence of the gas description 

We start by varying the gas parameters. This is discussed 
for the most general class of the DM, i.e. Model III (allowing 
for both convection and reacceleration) . For obvious reasons 
(see below), our conclusions also hold for Models 0, I, or II. 

In ID models, the surface density of the gas in the 
model corresponds to the average of the true gas surface 
density (which depends on th e position in the Galaxy ) over 
the effective diffusion volume (jTaillet fe Maurinll2003l ) . This 
is a reasonable assumption to make for stable nuclei (see 
Paper II). However, the details of the volume over which 
to calculate this average are not straightforward. Moreover, 
even if a more realistic distribution of gas were to be used, 
the latter is not free of uncertainties. We thus allow for 
some uncertainty in this input. We also vary the fraction of 
hydrogen relative to helium. 



4.1. Influence of the hydrogen fraction 

In Table [3J the second and third lines (compare with the 
first line) show the effect of changing the hydrogen frac- 
tion in the ISM: the transport parameters are changed at 
most by ~ 5%. There is a systematic trend for 6 to de- 
crease with smaller fractions of helium in the gas. However, 
the uncertainty on the hydrogen and helium fraction is not 
more than a few %. We therefore conclude that this has no 
strong impact on the derived transport parameters. 



4.2. Influence of the surface density 

The fourth and fifth lines show the effect of changing the 
surface gas density Sism = 2/migM- Whereas it has a small 
impact on the diffusion slope 5, the other transport param- 
eters arc strongly affected. The change in the parameters 
can be understood if we look at the grammage of the DM 
(IMaurin et al.ll2002|) . In t he purely diffusive regime, we have 



fe.g.. IMaurin et al.ll2006l) 



pure— D M 



Zi S MmvL 
IK 



(7) 



where fh is the mean mass of the ISM. Let Kq , S ro1 , V* , 
and V* el be the best-fit parameters obtained for a surface 
density of gas 57 ef (first line of the Table) . We remind that 
L is fixed here. If the surface gas density is rescaled, i.e. 
E new = x x E rof , in order to keep the same grammage in 
Eq. ([7]), we need to have K\ 
get in Table [2] 

The same reasoning holds for the convective wind. In 
presence of a consta nt wind, the full ex pression for the 
grammage reads (e.g.. IMaurin et al.|[2006h 



'new 
L 



x x K™ 1 . This is what we 



2V r . 



(8) 



With the above rcscaling for E 

x x 



:w and^ ow . 



V c . Again, this is recovered in Table [5J 



we get V c n 
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Table 2. Best-fit transport parameters for different ISM. 



Gas 

E 


/h 


A'o bcst x 10 2 
(kpc 2 Myr" 1 ) 


^bcst 


(km s- 1 ) 


u best 
' a 

(kms- 1 


X 2 /d.o.f 

) 


E ref 


rrcf 

/h 


0.48 


0.86 


18.8 


38.0 


147 


E ref 
E rof 


95% 

80% 


0.53 
0.41 


0.83 
0.90 


18.6 
19.3 


38.1 
37.7 


1.24 
2.14 


1 v-ircf 

2E ref 


/■rcf 

Jh 

/■rcf 


0.25 
0.92 


0.85 
0.86 


9.5 
37.3 


19.4 
74.6 


1.45 
1.51 



Note 2. Model III for L = 4 kpc. For all settings, we keep fixed 
(n e ~ ) = 0.033 and T e ~ 10 4 K. The reference parameters for the sur- 
face density S ref and the hydrogen fraction f^ 1 are given in Tabic [T] 

The last parameter is the Alfvenic speed, for which we 
need to consider the whole transport equation. The speed 
V a is used to define K pp [see Eq. ©], which appears in the 
terms b and c [Eqs. ([3]) and (TJJ] of Eq. ([1]). If we consider 
the latter equation and the transport operator C in Eq. © , 
the rescaling E now = x x S rof leads to £ ncw = x x £ ref . This 
implies 6 now = xx & rof , c ncw =n c rof , and S now = xx S Tct . 
Let us consider the three terms each in turn. For b, if we 
look at Eq. ([3]), this x factor is automatically ensured for 
the Coulomb and ionisation losses (they depend on the gas 
surface density), and also for adiabatic losses (as V^ cw = 
x x V c rcf ). This implies K™ w = xx K™ 1 . From Eq. ©, we 
have K pp oc V^/Ko, such that the previous equality yields 

(C W )7^ CW = x x {V^ f/Kf. 

This gives V a ncw = x x V a (the same information is given 
by the c term) , in agreement with the results of Table [5] 

4.3. Summary for the gas surface density 

An uncertainty of x% on the gas surface density leads to 
SystUnc of x% on Kq, V c and V a , but it leaves the diffusion 
slope 5 unchanged. This is true, whatever the value for the 
halo size L. For the sake of completeness, we should note 
that this also affects the source term: the source abundances 



5. Influence of the source spectrum 

In Paper I, we showed that some correlations exist between 
the source spectra and the transport parameters. Below, we 
investigate in more detail how the transport parameters are 
changed if we assume different input for a (spectral index 
of the sources) and r/s (related to the spectral shape at low 
energy: 

Q{E) oc ^ s TZ- a . (9) 
5.1. Influence of the source slope a 

It is usually said that the source spectrum energy depen- 
dence Q(E) factors out of the secondary-to-primary ratio. 
This agrees with the results gathered in Table [3J which 
shows the effect of varying the source slope a (or equiva- 
lently the slope of the measured fluxes 7 = a + 6). This 
is especially true for Model I (no rcacceleration) , whereas 
if reacceleration is included (Models II and III), a few % 



Table 3. Best-fit transport parameters for different slopes 
7- 



7 


K$ CBt x 10 2 
(kpc 2 Myr^ 1 ) 


^best 


tt best 

(kms^ 1 ) 


^best 

(km s" 1 ) 


x 2 /d.o.f 


0: 2.55 


4.13 


0.39 






29.7 


0: 2.65 


4. 08 


0.40 






28.8 


0: 2.75 


4.02 


0.40 






27.9 


0: 2.85 


3.97 


0.40 






27.0 


I: 2.55 


0.42 


0.93 


14.0 




11.4 


I: 2.65 


0.42 


0.93 


13.5 




11.2 


I: 2.75 


0.42 


0.93 


13.1 




11.0 


I: 2.85 


0.42 


0.93 


12.7 




10.9 


II: 2.55 


9.41 


0.24 




74.2 


4.52 


II: 2.65 


9.76 


0.23 




73.1 


4.73 


II: 2.75 


10.1 


0.23 




72.1 


4.94 


II: 2.85 


10.4 


0.22 




71.0 


5.15 


III: 2.55 


0.52 


0.84 


18.9 


39.7 


1.25 


///: 2.65 


0.48 


0.86 


18.8 


38.0 


1.47 


III: 2.75 


0.45 


0.87 


18.8 


36.4 


1.75 


III: 2.85 


0.42 


0.88 


18.6 


35.0 


2.10 



Note 3. Models 0, I, II, and III for L = 4 kpc. The parameter 7 is 
the slope of the propagated fluxes (7 = + 6). The source spectrum 
is Q(E) oc pvsn- a , and -q s = -1. 



change on the parameter 7 leads to a few % change on Va, 
-K01 an d S. 

From the last column of the same table, we see that 
on a statistical basis, the best model is Model III, then 
Model II, Model I, and finally Model 0. The relative merit 
(goodness-of-fit) of each class of models (III, II, I, and 0) 
remains unchanged when a is changed. 

5.2. Influence of the low-energy source spectrum through 775 

The second parameter that can be varied in the source term 
is the low-energy shape /3 ?,s [see Eq. ([5])]. which is a priori 
not very well constrained, both from observational and the- 
oretical points of view. The evolution of the best-fit trans- 
port parameters as a function of r/s is shown in Fig. [TJ for 
Model (dash-dotted lines/stars), I (dotted lines/plusses), 
II (dashed lines/crosses), and III (solid lines/open circles). 

The three configurations have marked minima in their 
X^i n (bottom panel), but not at the same r]s- Models 
and I have a similar x 2 dependence (weak impact of the 
presence of the wind), as have Models II and III (both af- 
fected by reacceleration). Understanding the exact depen- 
dence of the transport parameters with 775 is not the goal 
of the analysis. The important point to underline is that 
changing r)s does not allow reconciling convection Models I 
and III with smaller i5. As i5 varies only in the range 0.3 — 0.5 
for Model (pure diffusion), we conclude that for present 
B/C data, large 8 are associated to the constant Galactic 
wind. The largest variation is observed for Model II (reac- 
celeration only), for which 5 is allowed to be either smaller 
or larger than <5j b j cst = 0.33 (associated with rjg cst w 1) for 
a smaller or a larger 775. For any 775, Model III always per- 
forms better than the other two models. 
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10- 





Influence of r| (source) - 

' " - -x - - * - - x. _ _ K „ _ — = 




r , , 1 , , , 1 , 





0.5 



80 

-r 60 

a 40 

>" 

20 

1 10 



I_l 1 1 1 1 1 1 1 




- — e— Model III 

*™ .MpdeJ^n^o^nv^) 
~ Model I (no reac.) 
_ Model (pure^dtf£)._ K ---> 
1 1 .*--**!<*" 1 1 1 1 1 


- — x- _ 







— x-"*'" 




» I) / 1 1 1 






? g - 




Fig. 1. From top to bottom: best- fit values of <5, V'a, 
V c , and the associated x 2 /d.o.f. as a function of ijs [see 
Eq. ©]. Here, 7 = a + (5 = 2.65. The vertical grey-dotted 
line is a guideline for the results at the default configuration 
Vs = -1- 

Table 4. Best- fit values from a simultaneous fit on B/C 
and O. 



Model 


J.B/C+O) 

'Is 


a (B/C+0) 


S (B/C+0) 


Xmin/d.O.f. 





-0.82 


2.22 


0.41 


25.3 


I 


-1.23 


2.23 


0.98 


9.23 


II 


1.24 


2.24 


0.35 


3.62 


III 


1.56 


2.30 


0.95 


2.31 



Note 4. The data used arc IMP7-8, Voyager 1&2, ACE, HEAO-3, 
Spacelab, and CREAM04 for the B/C ratio, and HEAO-3 for the 
oxygen flux. 



5.3. Fitting a and rjs using a primary flux 

The source parameters a and rjs are not completely free. 
They can be fitted directly on primary elemental fluxes (sec 
Paper I) . If we determine simultaneously the transport and 
source parameters from a fit on B/C + O HEAO-3 data, we 
obtain the best-fit values gathered in Table SJ For Models II 
and III, Tig^ C+ ° is not far from rj ] g est found from the B/C 
fit (see Fig. [T]). For Model II, using the above constraint 

(Vs /C+ ° = 1-24) selects the diffusion slope 5 B / c +° = 0.35, 
close to the Kolmogorov spectrum for the turbulence. 



The surprising result is that the value for a B / c+ ° is 
quite resilient to any slope obtained for 6. This is odd, 
since it leaves the propagated slope 7 = a + 5 quite uncon- 
strained, leading to either 2.6 (Model II) or 3.2 (Model III). 
This indicates that in Model III, pure diffusive transport 
(i.e. 7 = a + S) has not yet been reached at TeV/n ener- 
gies. 

5.4. Summary for the source effect 

The dependence of the transport parameters with the 
source parameters a and rjs are very different. We checked 
that the transport parameters do not significantly depend 
on the source parameter a (in particular, S is unaffected). 
Still, correlations exist between the source and the other 
transport parameters (see Paper I) , and from Table O we 
can estimate that they affect their best-fit values at most 
by a few ten percent. 

On the other hand, Fig. [T] shows that the transport pa- 
rameters (including 5) are quite sensitive to the low-energy 
source spectrum parametcriscd by rjs- In any case, this does 
not allow reconciliation of Model III with small 5. Taking 
a sizeable range for rjs always lead to S > 0.6. This is even 
worse for Model I (no reacceleration). At variance, Model II 
(no convection) is very sensitive to 775. However, 775 is con- 
strained by primary fluxes measurements, and by no means 
could it span the range shown in Fig. [TJ When this con- 
straint is considered, the best-fit 5 for Model II turns out 
to be 0.35. Given the sensitivity of the latter to 77s and the 
possible bias in low-energy data (because of the solar mod- 
ulation) , this is consistent with a Kolmogorov spectrum. If 
the constraint of the primary flux is taken into account, the 
uncertainty on the transport parameters (associated to the 
source parameter r/s) can be estimated as ~ 10 — 20%. 

The puzzling point is that, although Model III is pre- 
ferred over Model II on a statistical basis, the latter seems 
to agree more with the expected theoretical value of S. This 
is more of further if we consider the quantity 7 = a + S, 
the asymptotic slope reached in the purely diffusive regime. 
The TRACER exp eriment finds 7 « 2.65 for all nu- 
clei (jAve et al.ll2008l ). whereas the bes t-fit slope from high - 
energy H and He data is 7 2.85 (|Donato et al.l l2009f) . 
From the propagated slope's point of view, Model II is 
again favoured over Model III, since the former leads to 
7 = 2.59, whereas the latter leads to 3.25 (read off TablcS]). 
If Model III is confirmed, as already emphasised above, this 
would mean that even high-energy data have not reached 
the purely diffusive transport yet. 

6. Influence of the low-energy diffusion coefficient 

We now turn to the effect of the low-energy shape of the 
diffusion coefficient. Assuming different diffusion schemes 
lead to different forms of the spat ial and mo mentum diffu- 
sion c oefficient (|Schlickeiserll2002h . Recently. iPtuskin et al.l 
(|2006l ) also argue that the form of the spatial diffusion coef- 
ficient can change at low energy, due to the possibility that 
the nonlinear MHD cascade sets the power-law spectrum of 
turbulence. 

We first consider several theoretically-motivated diffu- 
sion forms given in the literature. They mostly differ at low 
energy, so that it is useful to rewrite Eq. (J5J) as 



K{E) = ^ T ■ K Q K S 



(10) 
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Table 5. K(E) and K pp for different schemes. „ 



Type of turbulence 



i It 



4/3p 2 V a 2 



LBI 


Leaky Box Inspired 




i 


S (4- 


<5 2 ) (4-6) 


SA 


Slab Alfven 1 




1 


S (4- 


^)(4-5) 


IFM 


Isotropic fast magnetosonic 2 — 5 






Mix 


Mixture S A and IFM 1-5 


P 1 - 






0.5 



Note 5. The spatial diffusion coefficient is K*_ 



01T -Ko -n" 




• Model HI 
Model II (no conv.) 
Model I (no reac.) 
Model (pure dif^.) 



Table 6. Best-fit transport parameters based on different 
low-energy dependence of the diffusion coefficient. 



Type 



K% CBt x 10^ 
(kpc 2 Myr" 1 



V 



(kms" 1 ) (kms- 1 ) 



XVd.o.f 



0: LBI 
0: SA 
0: IFM 
0: Mix 



LBI 

SA 

IFM 



3.48 
4.08 
4.30 
3.71 

0.40 
0.42 
0.42 



0.45 
0.40 
0.38 
0.43 

0.94 
0.93 
0.93 



13.6 
13.5 
13.5 



17.5 
28.8 
36.7 
23.7 

12.0 
11.2 
11.6 




I: Mix 


0.41 


0.94 


13.5 




12.0 




t 1 1 1 1 | 1 1 


1 1 1 1 1 | 1 1 


1 1 1 1 1 1 i 
J,...-*-""?* 




II: LBI 


5.50 


0.38 




65.0 


1.61 


/d.o.f 




* H * 


I 6 "- **■■■■+■■■■ + ■■■■+-....+ + .... + ..„, 

; , ^ " 


.X* 


//: SA 
II: IFM 


9.76 
14.0 


0.23 
0.16 




73.1 
18.9 


4.73 
6.86 


1 


I", , T T 1 . 1 


1 1 \ T 1 1 1 1 T 1 


.X** k 


1 1 — 


II: Mix 


7.13 


0.32 




12.8 


2.03 




4 -3 


-2 -1 


1 


\ 



III: LBI 

///: 

III: IFM 
III: Mix 



0.70 
0.48 
0.49 
0.73 



0.78 
0.86 
0.85 
0.77 



18.0 
18.8 
18.9 
17.8 



47.1 
38.0 
45.6 
57.4 



0.87 
1.47 
1.25 
0.93 



Fig. 2. Same as Fig. [I] but now as a function of t\t [see 
Eq. (flU)) ]. The vertical grey-dotted line is for the default 
configuration t\t — 1 (SA). 



Note 6. Models 0, I, II, and III for L = 4 kpc. SA corresponds to the 
reference DM used throughout the paper. 

In the second step, we let t\t vary over a wide range in 
order to draw more general conclusions. 

6.1. Influence of the turbulence scheme 

A few turbulence schemes are gathered in Table EJ The 
associated best-fit values of the transport coefficients are 
presented in Table El When only convection is present 
(Model I), the low-energy form of K(E) is irrelevant, as 
seen in Table [BJ For models with reacceleration (Models II 
and III), it significantly affects almost all parameters. If 
there is no wind (Model II), the effect is maximum on 8. 
The model with pure diffusion (Model 0), or with both con- 
vection and reacceleration (Model III) , falls in-between. For 
Model III, the cases SA and IFM, on the one hand, and LBI 
and Mix, on the other, give very similar results. This is eas- 
ily understood as the quantities 7j|A = 2 — 5 and 77^ F = 1 
(respectively T]^ B1 = 1 — S and 77^ lx = 0) are roughly equal 
for 5]^' ~ 1. 

The important result is that Model III is mildly sensitive 
to the diffusion scheme, with an uncertainty of a few ten 
percent scatter on all the transport parameters. However, 
Model II parameters are extremely sensitive to the diffusion 
scheme (more than a factor of 2 scatter). Depending on 



the case considered, 5 is found in the range 0.16 — 0.38. 
The hierarchy of Xmin among the various models is always 
conserved. 

6.2. Generalisation to any ifr 

We generalise the analysis by allowing for any value of the 
parameter r/T in the diffusion coefficient K(E). In doing so, 
we do not seek to provide sound physical motivations for 
the range tested. In this section, whatever the value of t/t, 
the diffusion coefficient in momentum space is assumed to 
follow K PP K KX = (4/3)p 2 V*/(5 (4 - S 2 ) (4 - 6)). 

The best-fit transport parameters and Xmin evolution as 
a function of rjx are plotted in Fig. [21 For tjt < —2, all the 
models converge slowly towards purely diffusive models (no 
convection, no reacceleration). But this is at the cost of a 
bad x 2 ( see bottom panel) . Based on the x 2 criterium, high 
values of 77T (> 2) are also disfavoured. The four configu- 
rations have marked minima in their Xmim corresponding 
to ?7^ cst w -2.75, -2.5, -0.25, +0.25 for Models 0, I, II, 
and III respectively, for which S bcst w 0.6, 0.6, 0.4, 0.8. For 
rj T —2.5, all models point to 6 ~ 0.5, close to a Kraichnan 
spectrum for turbulence. 

For Model III, for a well-chosen value of -t/t, the diffusion 
slope could be decreased at most down to 5 = 0.5, but such 
a configuration does not correspond to the minimal x 2 f° r 
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this model so is excluded. For Model II, almost any value 
of S can be reached. 

6.3. Summary for the diffusion coefficient effect 

Similars to the effect of the source parameter r/s, varying 
the parameter t]t on a wide range i) does not to allow 
S < 1/2 to reached for Model III, and ii) strongly affects S 
for Model II. For the latter, the best-fit t\t leads to 5 *=s 0.4, 
slightly more than the Kolmogorov spectrum of turbulence. 

From a statistical point of view, Model III is again pre- 
ferred, except in the region — 3 < t?t < —1.5 where V c drops 
to zero, so that it is equivalent to Model II. For t\t < —2.75, 
V a also drops to zero, so that all models are close to the 
purely diffusive case (Model 0), which favours 6 ~ 1/2. As 
before, one of the interesting features of Model II (reac- 
celeration only) is its versatility. As the range of allowed 
T]t remains unspecified to some extent, Model II does not 
point to any specific value of i5. It can accommodate values 
as low as 1/3, going through 1/2 and even higher values, if a 
sharp turn-off ( negative values of rfr ) exists in the diffusion 
coefficient (e.g.. iPtuskin et aDl2006h . 



Table 7. Best-fit transport parameters for various cross- 
section sets. 



X-files A' bcst x 10 2 V^ 1 x'/d.o.f 



(kpc 2 Myr" 1 ) 




(km s _1 ) 


(km s _1 ) 




0: WKS98 


2.66 


0.53 






23.5 


0: SOI 


3.40 


0.40 






34.5 


0: W03 


4. 08 


0.40 






28.8 


0: GAL09 


3.83 


0.46 






28.4 


I: WKS98 


0.45 


0.95 


10.4 




12.0 


I: S01 


0.25 


1.01 


11.9 




13.5 


I: W03 


0.42 


0.93 


13.5 




11.2 


I: GAL09 


0.49 


0.95 


13.6 




12.3 


II: WKS98 


7.19 


0.31 




71.5 


3.40 


II: S01 


9.27 


0.22 




68.8 


6.56 


II: W03 


9.76 


0.23 




73.1 


4.73 


II: GAL09 


10.0 


0.26 




85.0 


4.03 


III: WKS98 


0.69 


0.80 


15.8 


43.5 


1.11 


III: S01 


0.25 


0.98 


16.2 


27.9 


3.01 


III: W03 


0.48 


0.86 


18.8 


38.0 


1.47 


III: GAL09 


0.65 


0.82 


21.7 


49.4 


1.53 



7. Influence of the production cross-sections 

Many reaction channels are required to calculate the pro- 
duction of secondary species. Semi-analytical formulae, 
semi-empirical approaches, or even fit to the data are cur- 
rently used in the literature to obtain the full set of pro- 
duction cross-section in CR physics. We show first the dif- 
ferent best-fit values obtained for different sets available in 
the literature. We then inspect what would be the effect 
of energy-biased cross-sections on the transport parameter 
derivation. 



7.1. Using different sets of fragmentation cross-sections 

We use four different sets of cro ss-sections. Among t hem, 
the most up-to-date are W03 (jWebber et al.l 120031) and 
GAL09 (from the GALPROP codeQ). The former is based 
on the semi-empirical approach of Webbe r and coworkers, 
initiated in the 90 's (jWebber et al.lll990D . The latter set 
takes advantage of the former approach, rcnormalising some 
cross-sections to sel ected data. The two others, WKS98 
(|Webber et al.l Il998t ) and S01 (private communication of 
Aime Soutoul, 2001), are a bit outdated, but they serve as 
an illustration here. 

The effect of these different sets on the parameter de- 
termination is shown in Table [7] The scatter in the values 
of the different parameters is not the same depending on 
the model considered. WKS98 requires less convection than 
any other set, and GAL09 needs more reacceleration than 
the others. Otherwise, it is difficult to find clear trends. For 
instance, the W03 and GAL09 sets give very similar results 
for Model I, but not as similar as for Models II and III. 

If we discard the S01 set (which is based on an un- 
published preliminary analysis), we can conclude that the 
SystUnc on the transport parameters related to the choice 
of the production cross-section is < 20%. This figure is sim- 
ilar to the uncertainty quoted for the parameterisation of 
the cross-section themselves (~ 10 — 20%). 



Note 7. Models 0, I, II, and III for L = 4 kpc. The four base sets of 
cross-sections are WKS98 from Webber et al. (1998), S01 from a pri- 
vate communication of Aime Soutoul (2001), W03 from lWebber et al.l 
<2003h . and GAL09 from the GALPROP code v50.1p (2009). 



7.2. Influence of a systematic energy bias in W03 
cross-sections 

A striking re sult of the updated cross-section formula- 
tion of W03 ([Webber et alj I2003D compared to WKS98 
(jWebber et al.lll998l) was a systematic energy bias (smaller 



1 http://galprop.stanford.edu/wcb_galprop/galprop_home.html 



cross-sections at higher energy). To address the effect of a 
possible residual bias, we allow for a systematic energy bias 
in all production cross-sections, either at low energy (LE) 
or at high energy (HE). Biased sets are obtained from the 
reference set W03 by applying a factor (E k / n /10)~ x below 
10 GeV/n (LE bias), and (E k/n ) x above 1 GeV/n (HE bias, 
but left constant above 10 GeV/n). This is illustrated for 
the 12 C— )• ( 10 B+ 11 B) cross-section, as shown in Fig. [21 As 
seen from this figure and the experimental data, biases x 
as large as 0.05 are marginally consistent with the data at 
low energy, but are acceptable at high energy. 

In Table El we only report the results of the high-energy 
biases on the best-fit values of the transport parameters. We 
checked that the result of a low-energy or a high-energy bias 
gives the same trends. Quite naturally, there is a significant 
correlation of the bias with the transport parameters. We 
note that if 5 increases, all the other transport parameters 
Kq, V a , and V c decrease. In principle, we would expect that 
5 increases if the bias is positive: larger cross-sections at 
high-energy produce more secondaries at high-energy, re- 
quiring a larger S to match the same B/C data. This is 
observed for Models and II. However, the reverse effect is 
observed for Models I and III, i.e., for models with connec- 
tive wind. 

Such an exercise has limitations since it is not realistic 
to expect such systematic energy biases for all production 
channels. However, it gives some trend and further confirms 
that Model III cannot be reconcile with small 5. 
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Fig. 4. Left panel: Model II (pure reacceleration). Right panel: Model III (reaccclcration and convection). PDF of the 
transport parameters (for the reference inputs) as obtained in Paper II, along with the best-fit values for different 
production cross-section sets. 




Fig. 3. Prod uction cross-sectio n for 12 C+H— i> 10 ' n B 
(adapted from I Webber et aDl2003h . The standard sets arc 
shown as solid lines (WKS98: red dots; GAL09: red down 
triangles; W03: black stars), and the biased sets in dotted 
(\x\ = 0.02) and dashed (|x| = 0.05) lines. 



7.3. Summary for the influence of fragmentation 
cross-sections 

The results when using different production cross-section 
sets are difficult to interpret. To check that the observed 
trends do not come from an error in t he code, we compare d 
our results with the LBM analysis of lWebber et al.l (|2003f) . 
These authors find a significantly smaller escape length de- 
pendence on rigidity (from P 6 to P 5 ^, going from the 
WKS98 to the W03 set. From both an LBM analysis with 
our propagation code and a DM analysis using the LBI case 
(see Table O with the addition of a rigidity cut-off (not dis- 
cussed in this paper), we also find this decrease of AS = 0.1. 



Table 8. High-energy biases on the W03 cross-section set. 



W03 

HE bias x 


^ best x 1Q 2 

(kpc 2 Myr" 1 


^best 


t^/ best 

(km s~ 


■^/-best 

J ) (kmV 1 


x 7d.o.f 

) 


0: 


+0.05 


4.05 


0.45 






26.5 


0: 


+0.02 


4.06 


0.42 






27.7 


0: 


+ 0.00 


4.08 


0.40 






28.8 


0: 


-0.02 


4.09 


0.38 






30.1 


0: 


-0.05 


4.13 


0.34 






32.5 


I: 


+0.05 


0.61 


0.89 


13.7 




12.1 


I: 


+0.02 


0.49 


0.91 


13.6 




11.4 


1: 


+ 0.00 


0.42 


0.93 


13.5 




11.2 


I: 


-0.02 


0.35 


0.95 


13.4 




11.2 


I: 


-0.05 


0.27 


0.98 


13.1 




11.8 


II 


+0.05 


9.20 


0.28 




78.0 


3.26 


II 


+0.02 


9.51 


0.25 




75.0 


4.13 


// 


■ +0.00 


9.76 


0.23 




73.1 


4.73 


II 


-0.02 


10.0 


0.21 




71.2 


5.35 


II 


-0.05 


10.4 


0.19 




68.1 


6.27 


II 


[: +0.05 


0.98 


0.74 


19.0 


50.0 


1.01 


II 


[: +0.02 


0.75 


0.81 


19.0 


42.7 


1.22 


// 


/: +0.00 


0.48 


0.86 


18.8 


38.0 


1.47 


II 


I. -0.02 


0.35 


0.91 


18.5 


33.4 


1.88 


II 


t: -0.05 


0.21 


1.01 


17.8 


26.6 


2.94 


Note 8. Models 0, I, II, 


and III 


for L 


= 4 kpc. 


The vari- 



ous sets are obtained after multiplying the W03 set by (Ek/n)^ 
above 1 GeV/n (and kept constant above 10 GeV/n), with x = 
+0.05, +0.02, 0, -0.02, -0.05. (corresponding to bottom to top 
curves on the right-hand side of Fig- [3jl . 

This is also confirmed for the more standard Model (pure 
diffusion) where we go from 0.53 to 0.40. 

We thus have to conclude, as seen from Table [7J that 
the effect of these cross-section sets is not the same for dif- 
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Table 9. Influence of biasing HEAO-3 data. 




IMP7-8 (0=250 MV) - Garcia-Munoz et aL 1987 
Voyagerl&2 ($=225 MV) - Lukasiak 1999 
ACE 97-98 (#=225 MV) - De Nolfo et al. 2006 - 
HEAO-3 (#=250 MV) - Engelraann et al. 1990 ' 
Spacelab (#=200 MV) - Swordy et al. 1990 
CREAM 04 (#=425 MV) - Ahn et al. 2008 



0.1 



0.05- 



Biased HEAO-3 data 

— x = +0.10 

— x= +0,05 
- x = -0.05 

i I 

1 10 



ElAffi (GeV/n) 



Fig. 5. Demodulated B/C data used in this paper (de- 
noted dataset F). Biased HEAO-3 data are shown as solid, 
dashed, and dotted lines. 



ferent classes of models: it strongly depends on the presence 
or absence of convection and, to a lesser extent, on reac- 
celeration. Nevertheless, for all models, the typical scatter 
in the best- fit values is a factor of 2 for Kq, ~ 50% for V a , 
~ 10% for 5, and ~ 5% for V c . Figure 0] illustrates that un- 
certainties in some input parameters (here the production 
cross-sections) provide larger SystUnc on the transport pa- 
rameters than their StatUnc calculated using the MCMC 
technique (taken from Paper II). 

As for the other effects inspected in previous sections, 
Model III always provide 5 larger than ~ 0.6. Model II 
provides values of S in the correct range for a Kolmogorov 
spectrum, but it never does better than Model III in terms 
of the goodness-of-fit. 



8. Influence of an energy bias in HEAO-3 B/C data 

It is obvious that most of the results derived in this pa- 
per rely on the confidence we place in the HEAO-3 data, 
since they are the most constraining data to date. Their 
very small error bars at high energy could over-constrain 
6. Several other studies use large error bars for high en- 
ergy B/C HEAO-3 points, c ontradicting the prescr iption 
given in the original paper (jEngelmann et all 11990ft . and 
this could bias their results. As an exercise, we allowed for 
an energy bias in these data. This approach is disputable so 
we do not wish to defend it strongly. We inspected whether 
a small bias in the data strongly affects the best-fit values 
of the transport parameters. 

The HEAO-3 data are biased above E c , using 
(E^/ n /E c ) x for £? k / n > E Cl with E c = 2.9 GeV/n, as shown 
in Fig. [5] As in previous sections, we fit Models 0, 1, II, and 
III on dataset F (as shown in Fig. [5]), where the HEAO- 
3 data have been replaced by the biased ones. The results 
have been gathered in Tabled Intuitively, in a similar fash- 
ion as for the bias in the cross-sections, we expect that a 
higher value of the B/C data at high energy (positive bias) 
leads to a smaller S. This is what is observed for all models. 
The effect of the bias on S is at its maximum for the pure 
diffusion model (Model 0), where AS ~ 0.2. For Models II 
and III, the effect is less pronounced owing to the pres- 



HEAO-3 
HE bias x 


Ko CBt x 10 2 
(kp^Myr -1 ) 


^best 


T7 bCSt 

' c 

(km s~ 


^best 

x ) (kmV 1 ) 


y 2 /d f 


0: +0.10 


4.97 


0.28 






22.5 


0: +0.05 


4.50 


0.34 






25.2 


0: +0.00 


4.08 


0.40 






28.8 


0: -0.05 


3.69 


0.46 






33.0 


t- -Hi 1 n 




n si 

U.Ol 


10.1 




Q 77 


I: +0.05 


0.44 


0.87 


13.5 




10.2 


I: +0.00 


0.42 


0.93 


13.5 




11.2 


I: -0.05 


0.40 


0.99 


13.6 




12.5 


II: +0.10 


10.1 


0.17 




56.6 


4.54 


II: +0.05 


9.95 


0.20 




64.8 


4.27 


II: +0.00 


9.76 


0.23 




73.1 


4.73 


II: -0.05 


9.56 


0.27 




81.5 


5.79 


III: +0.10 


0.31 


0.85 


19.1 


30.0 


1.32 


III: +0.05 


0.42 


0.83 


19.0 


34.8 


1.16 


III: +0.00 


0.48 


0.86 


18.8 


38.0 


1.47 


III: -0.05 


0.50 


0.90 


18.7 


40.0 


2.18 



Note 9. Models 0, I, II, and III for L = 4 kpc, where HEAO-3 
data are biased using the formula (Ek/ n / E c ) x for E k / n > E c , with 
E c = 2.9 GeV/n and x = +0.1, +0.05, 0., -0.05. 

ence of reacceleration, with a change of only A<5 ~ 0.1. 
However, Model III (convection and reacceleration), which 
is statistically preferred over Model II, still points to un- 
comfortably high values for S. The only way out would be 
to have a non-systematic energy effect in the data. For in- 
stance, t he result on the fit to the recently published AMS- 
01 data (|Tomassetti fc AMS-01 Collaboratio"n1l2009ft leads 
to S Ki0.5 (see Paper II). 

Finally, for completeness, we checked that the effect of 
the solar modulation parameters (we took A(f> ± 50 MV for 
the HEAO-3 data) was negligible on all the derived trans- 
port parameters. 



9. Summary and discussion 

The guiding questions for this paper were the following. 
First, do uncertainties in the input ingredients lead to 
systematic uncertainties (SystUnc) on the derived trans- 
port parameters greater than their statistical uncertainties 
(StatUnc)! Second, can we reproduce the various values of 5 
given in the literature, and which one should be preferred? 
We elaborate on the answers below before concluding. 



9.1. SystUnc and StatUnc: which ones dominate? 

From the analysis of Paper II, the StatUnc found for the 
various transport parameters typically fall in the range 
5% — 10% (see Fig. @|. The SystUnc generated by each 
input we varied are summarised below. 

— A change of x% in the gas surface density Eism trans- 
lates into a change of x% for all transport parameters 
but <5, which is unaffected. Such a simple scaling can be 
used for comparing models in the literature, as different 
Sism are generally used. Depending on how confident 
we are in the measurement of Sism, we may conclude 
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that, for this input, the generated SystUnc is similar to 
or larger than the StatUnc. 

— The transport parameters arc very sensitive to the low- 
energy spectral shape of the source parameter rjs, but 
not to the source slope a. Both parameters can be de- 
termined by including data of a primary flux in the fit. 
When doing so, the SystUnc and StatUnc are of the 
same order. 

— If we let vt — which paramctcrises the low-energy 
shape of the diffusion coefficient — free in the range 
[—2., 1.], the SystUnc completely dominate the StatUnc. 
Diffusion coefficients with a sharp t urn-off at low-energy 
are possible, as shown in Fig. 1 of iPtuskin et al.l 120061 
(their dashed line corresponds to our t\t = —2). 

— For the cross sections, as seen in Fig. 2J the SystUnc 
can also be larger than the StatUnc for some of the 
transport coefficients. The use of other secondary-to- 
primary ratios, such as Li/C and Be/C, could help to 
cross-check the consistency of each cross-section set, and 
thus decrease the SystUnc. 

These conclusions hold for all classes of models (0, I, II, 
and III). 

9.2. Which value of 5 should we trust? 

The two key ingredients for detcrminaning 6 are i) the low- 
energy spectral shape of the diffusion coefficient as parame- 
terised by i]t [sec Eq. (fTQ]) ]. and ii) the presence or absence 
of a wind. These two effects allow 5 to reach val ues as low as 
0.2 up to 0.9. Such a wide range is also found in I Jones et al.l 
(|200lh . depending on the model they consider (convection, 
turbulent diffusion, stochastic reacceleration) . 

We recover S « 0.3 in the pure reacceleration 
case (no convection), as in other propag ation codes 
(jStrong fc Moskalenkolll998t I Jones et al.ll2001l i. As soon as 
a convective wind is included, the diffu sion slope is large 
(5 ~ 0.8), as obtai ned in other studies (| Jones et al.ll200lt 
iMaurin et aTll2001l ). The effect of the low -energy shape of 
the diffusion coefficient is illustrated by the fact that stan- 
dard LBM and standard diffusion models lead to different S 
(more details are given in Putze, Derome & Maurin, ICRC 
2009). For LBM, which are equivalent to di ffusion models 
if K( E) = K TZ 5 , we recover 6 « 0.5 - 0.6 ([Webber et all 
120031 Paper I. 

It is reassuring to see that the results of the various 
propagation codes used in the literature are consistent. 
However, this does not settle the question of which value of 
5 we should trust. In principle, one advantage of our anal- 
ysis is that it should clearly point to a preferred value for 
<5, by a mere comparison of the best x 2 values. Indeed, in 
almost all settings considered, Model III (convection and 
reacceleration) performs better than any other, the second 
best being Model II (no convection). This behaviour is il- 
lustrated by the Xmin dependence as a function of 5, as 
shown in Fig. [6] For a standard diffusion coefficient (black- 
solid curve), there is a clear minimum at high 5, and a 
less marked local minimum at low 5. This second minimum 
matches the Xmin °f Model II (black-dashed curve). The 
value of the Xmin slightly decreases when t\t decreases (red 
curves), then increases for t\t < (grey and violet curves). 
In the process, the best-fit value for 8 is displaced to smaller 
values, the local minimum at low S overtaking and merging 
with the high 5 one. 



6 

10 



1 



Fig. 6. Evolution of Xmin/d-O-f- as a function of S for III 
(thick lines with symbols) and II (thin lines). The four cases 
correspond to different r\T (low-energy shape of the diffu- 
sion coefficient). 



Table 10. Best-fit parameters for a few selected configura- 
tions. 



Model 


Vt 


K x 10" 


5 


V e 


V a 


X Vd.o.f 






(kpc 2 Myr" 




(krns" 1 ) 


(km s~ 




nit 


SA 


0.481 


0.856 


18.84 


37.98 


1.47 


m/ir 


-1.3 


3.161 


0.512 


0. 


45.35 


2.26 


1/0* 


-2.61 


2.054 


0.613 


0. 




3.29 




SA 


9.753 


0.234 




73.14 


4.73 



f Best-fit transport parameters for standard Model II and III. 
\ Best-fit parameters with t\t free (no reacceleration). 
* Best-fit parameters for fixed r\^. 



Note 10. Standard models refer to SA diffusion coefficients (see 
Table [5]). Alternative models show their different values for nr, 
as parameterised from Eq. (|10p . 

For standard diffusion schemes (SA), the best class of 
models is Model III. The main problem is that whatever 
ingredients are changed in this model, we are always left 
with 8 > 0.6, and most of the time with even higher values, 
<~ 0.9. Such high values, always associated to the presence 
of the constant convective wind, are difficult to explain. 
Should we choose to discard this Model III, we would be 
left with the task of explaining why the statistical analysis 
fails in the context of CR physics. A way out is provided, 
as seen in Fig. [6j if we let free t\t (low-energy shape of the 
diffusion coefficient). Indeed, the only case where Model II 
(reacceleration only) performs as Model III is for t)t in the 
range —3 < t]t < —1.5. In such a configuration, the best- 
fit value for the convective velocity of Model III drops to 
zero, so that the best model only has reacceleration (i.e. 
equivalent to Model II): we term this class of models III/II. 
The next best-fit class of models is obtained if we discard 
reacceleration (Model I-like) and again let t\t go free. As 
for the previous case, the best-fit value for the convective 
velocity drops to zero: we denote this class 1/0. In terms 
of the best Xmim ^he standard diffusion Model II has only 
the fourth-best x 2 value. The best-fit values for these four 
classes of models are gathered in Table [TU1 sorted according 
to their figure of merit. 
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Fig. 7. Best- fit B/C ratio for standard Model II (thick red) 
and III (thick black) using the reference setting (solid and 
thin dashed lines are for the cross-section sets W03 and 
GAL09). The special best-fit Models O/I and arc 
also plotted in thin grey lines (light and dark shade respec- 
tively). 



The best-fit curves for these four relevant configura- 
tions are shown in Fig. [JJ From this graphic view, we can 
see that Model II matches the data best at very-low en- 
ergy, Model III at intermediate energies, and Models O/I 
and III/II at high energy. While we found that a system- 
atic bias in HEAO-3 data would not change the ordering 
of which is the best model (see Sect. [5]), we see from Fig. [JJ 
that a more complicated pattern in the data may be more 
effective in changing which model is best. 

To summarise, when a standard scheme for the diffu- 
sion coefficient is assumed, Model III is always preferred 
over Model II. However, the latter seems to agree more 
with theoretical expectations for 5. Model II is also more 
consistent with the measured propagated slope 7 = a + 6. 
Because of this oddity, we may not be able to give a defi- 
nite answer about the value of 5. This is further complicated 
when ?7t is left free; in that case, best models are always 
without convection, and even the pure diffusive model is 
redeemed (although both have difficulty reproducing the 
B/C peak at GcV/n energies). We see that many uncer- 
tainties show up at GeV/n energies. As also illustrated by 
the various predictions for various 5 in Fig[7J the higher the 
energy, the closer we can expect to reach the purely diffu- 
sive regime. As a result, high-energy B/C data are desired 
to unambiguously pinpoint the value of 6. 



value of 8 is real (in that case it needs to be explained theo- 
retically), or that any model with convection should be ex- 
cluded (which contradicts the fact that winds are observed 
in many galaxies) . Even if we adopt the last alternative, no 
firm conclusions can be drawn on the value of 5. Indeed, if 
the statistical analysis is relaxed, a large category of mod- 
els are redeemed, attaining any value for S between 0.3 and 
0.9. These models may be purely diffusive, with convection 
and/or reacceleration, and are very sensitive to the shape of 
the low-energy diffusion coefficient (which is not prescribed 
theoretically for the moment). Data at higher energy are 
needed to solve this question. More constraints can also be 
obtai ned by combining several secondary-to-primary ratios 
(e.g.. IWebbe^ll997allbt) . This is left for a later study. 

This study has limitations. For instance, we only varied 
the source and diffusion parameters according to simple pa- 
ramcterisations. More complicated dependences could have 
been inspected. However, it is worth recalling that it may be 
dangerous to introduce too many ad hoc prescriptions, be- 
cause the statistical meaning — already unclear when com- 
paring the different classes of models — becomes less and less 
obvious as the number of parameters and models tested 
increase. In the framework of homogeneous and isotropic 
diffusion coefficients, a maybe more important issue is the 
question of the Galactic wind. A constant wind was cho- 
sen because of the simplicity of the solutions (of the cor- 
responding diffusion equation). On the one hand, Galactic 
winds are ubiquitous. On the other, as shown in this study, 
constant wind cannot accommodate a realistic slope of the 
diffusion coefficient. A linear wind may provide different re- 
sults. We are implementing a numerical solver in our prop- 
agation code to inspect this issue. Another possibility that 
cannot be ruled out is that the HEAO-3 data suffer from 
non-trivial systematics. 

Finally, the starting point of the paper was a comparison 
between systematic uncertainties (generated by uncertain- 
ties in the input ingredients) and statistical uncertainties 
for the values of the transport parameters. That the for- 
mer can be larger than the latter shows that many efforts 
are still needed in CR physics, especially for the production 
cross-sections, before one can take full advantage of any sta- 
tistical analysis, as performed in Papers I and II with an 
MCMC technique. Some of these issues may be resolved as 
new data on cosmic-ray nuclei are being released (cream, 
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9.3. Conclusion 

In the past years, we have promoted and used a 
model favouring both co nvection and reacceleration (e.g., 
iMaurin et al.|[200ll l2002t and subsequent studies) from sta- 
tistical criteria. The main and known problem of this model 
lies in its uncomfortably high value for 8 (S ~ 0.8). Such a 
model is also preferred in the present study. In addition, we 
found that the high value for 6 is extremely resilient to any 
change in the setting, which leaves us with several alterna- 
tives: assume that there are complicated biases in the data 
that conspire to give high S in such models, that this high 
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